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ABSTRACT 

Low-mass binary millisecond pulsars (LMBPs) are born with very small 
orbital eccentricities, typically of order ~ In globular clusters, 

however, higher eccentricities ^ can be induced by dynamical interactions 
with passing stars. Here we show that the cross section for this process is much 
larger than previously estimated. This is because, even for initially circular 
binaries, the induced eccentricity ef for an encounter with pericenter separation 
Up beyond a few times the binary semi-major axis a declines only as a power-law, 
ef (x (rp/a)“®/^, and not as an exponential. We find that all currently known 
LMBPs in clusters were probably affected by interactions, with their current 
eccentricities typically an order of magnitude or more greater than at birth. 
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1. Introduction 


Low-mass binary millisecond pulsars (hereafter LMBPs) are thought to be formed 
when an old neutron star in a binary system accretes material from a red-giant companion. 
This leads to the production of a spun-up (recycled) pulsar with a low-mass white-dwarf 
companion in a wide, circular orbit. This formation process is well understood theoretically 
(for recent reviews, see Verbunt 1993; Phinney & Kulkarni 1994). The predicted relations 
between white-dwarf mass and orbital period (Joss, Rappaport, & Lewis 1987) and 
between orbital eccentricity and period (Phinney 1992) are in remarkable agreement with 
observations. The orbits tend to be very efficiently circularized by tidal torques and viscous 
dissipation in the red-giant envelope during the mass-transfer phase. Residual eccentricities 
Ci ~ 10“®-10“^ for LMBPs at birth are explained as arising from small fluctuations of 
the external gravitational held of the red giant caused by the convective motions in its 
envelope (Phinney 1992). In globular clusters, the wide orbits of LMBPs can be perturbed 
significantly by passing stars. Thus final eccentricities Cf ^ Ci could be observed today, 
and the measured values contain important information about the dynamical history of the 
binaries and their environment. In dense clusters binaries can even be completely disrupted 
through interactions, so we expect to And LMBPs with eccentricities in the entire range 
10 “® ^ ^ 1 . 

Using numerical scattering experiments, Rappaport, Putney, & Verbunt (1989) 
computed cross sections for inducing eccentricity in binaries with initially circular orbits. 
More recently, Phinney (1992) discussed the problem on the basis of approximate analytic 
expressions for the induced eccentricity derived from analytic perturbation theory. These 
two studies agree in predicting that the induced eccentricity decreases exponentially with 
increasing pericenter distance Up. In contrast. Hut & Paczyhski (1984), also on the basis 
of numerical work, found that the induced eccentricity decreased as a power-law, with 
Of (X and speculated that this was coming from a non-vanishing higher-order secular 
perturbation effect (At lowest order secular perturbation theory gives a vanishing result for 
the induced eccentricity when the binary orbit is initially circular). 

The purpose of this Letter is to shed new light on this question, by presenting and 
discussing our own analytic results for the induced eccentricity (§2). On the basis of detailed 
comparisons with numerical integrations, we believe these analytic results to be correct 
over the entire range of interest for LMBPs, 10“® ^ ^ 1. The previously published 

expressions in which ef decreases exponentially with Up can grossly underestimate the effect 
when ^ 0.01. Indeed, we And that over the intermediate range where ^ 10“^, the 

induced eccentricity generally varies as a power-law oc making the cross section 

much larger than would be estimated on the basis of previous work. The implications of 
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our results for LMBPs in globular clusters are discussed in §3. 


2. Eccentricity Perturbations 


For binaries with initial eccentricity ^ 0, it is easy to show, in lowest order, that 
secular perturbation theory gives 

= -”PW-cosU). (1) 

Here mi and m 2 are the masses of the binary components, m 3 is the mass of the passing 
star, Mi 2 = mi + m 2 , M 123 = M 12 + m 3 , a is the semi-major axis of the binary, H the 
longitude of the node (measured in the orbital plane of the binary from its pericenter), and 
i is the inclination (f = 0 is coplanar prograde; f = vr is coplanar retrograde). Throughout 
this paper we assume for simplicity that m 3 is initially on a parabolic orbit with respect 
to the center of mass of the binary. A first immediate consequence of equation (1) is that, 
except for extreme mass ratios, one always has |Ae|/ei <C 1 , as long as Tp is at least a few 
times a. Thus, within the range of applicability of equation (1) (which includes the Tp ^ 00 
limit; see below), encounters with passing stars produce a negligible change in eccentricity. 
Second, we see immediately that this lowest-order expression goes to zero for = 0. We 
must therefore ask if higher-order terms could not become dominant over some intermediate 
range of values for Vp. 


Indeed, a secular perturbation calculation at the next order, for = 0, gives 
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Here uj is the longitude of the pericenter of the third body, measured from the node. Notice 
that expression (2) vanishes in the equal-mass case (mi = m 2 ). A detailed derivation of 
expressions (2) and (3) will be given in Heggie & Rasio (1995; hereafter HR). The easiest 
derivation that we have found exploits the Lenz vector (or ’’eccentric axis”; cf. Pollard 
1976). If only the quadrupole interaction is included the result is equation (1), whereas the 
octupole interaction gives equations (2)-(3) when = 0. 


Equations (l)-(3) were obtained from secular perturbation theory, i.e., by averaging 
over the orbital motion of the binary. This is a valid assumption for sufficiently large Tp, 
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when the encounter is quasi-adiabatic. At smaller Up, however, the duration of the encounter 
can be comparable to the orbital period of the binary, and departures from adiabaticity can 
be important. Small deviations from adiabaticity can be calculated analytically using the 
method developed by Heggie (1975, §5.4) for estimating the changes in the semi-major axis 
of the binary (which vanish to all orders in secular perturbation theory). Using this method 
we obtain, for = 0, 
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with 
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( 5 ) 

Here (j) is the orbital phase of the binary at pericenter. Notice that this expression is 
maximum for a coplanar prograde orbit [i = 0), and vanishes for a coplanar retrograde 
orbit [i = vr). 


Our theoretical expressions (l)-(5) are in excellent quantitative agreement with the 
results of numerical integrations. Extensive comparisons will be presented in HR. Here, in 
Figure 1, we only illustrate this agreement for a typical LMBP case, and for f = H = a; = 0 
The exponential expression (4) becomes important for Cf 0.01, but the power-law 
expression (2) dominates for ^ ^ 0.01. The dependence on phase is quite strong in 

the exponential regime, but becomes completely negligible in the power-law region. When 
combining expressions (l)-(5), we find that the best agreement with numerical results is 
obtained when we simply use the expression giving the largest ef for fixed Up. 


3. Implications for Binary Pulsars 


An immediate implication of our results, evident in Figure 1, is that cross sections for 
induced eccentricity can be over an order of magnitude larger than would be estimated on 
the basis of the exponential result (eq. [4]). This is for parabolic orbits, or equivalently, 
hard binaries (for which the cross section a oc Up). In general, for slightly hyperbolic orbits, 
the induced eccentricity at fixed Up increases with relative velocity, and the relative increase 
in the cross section is even larger (Detailed analytic and numerical results for hyperbolic 
orbits will be presented in HR). 

We can easily estimate the timescale t>e for inducing an eccentricity Cf > e. For 
ei < e ^ 0.01, this is determined by setting the angle-averaged expression (2) equal to e. In 
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the parabolic limit, and averaging over a Maxwellian distribution of relative velocities, we 
find 

f>e ~ 4 X 10“ yi' (ei < e ^ 0.01) ( 6 ) 

Here 724 is the number density of stars (assumed for simplicity to be all of mass m 3 ) in units 
of 10 ^pc“^, Uio is the one-dimensional velocity dispersion in units of 10 kms“^, Pd is the 
LMBP ’s orbital period in days, and we have used the same typical values of the masses 
as in Figure 1 (the dependence on the masses in general is very weak over the relevant 
range for LMBPs in globular clusters). Deviations from the parabolic limit remain small for 
binaries with P ^ 10^ days. For e 0.01 a similar procedure using expression (4) gives 

t>e ~ 2 X 10“n4^nioPd~^^^ yr (e 0.01) (7) 


Equations ( 6 ) and (7) can be used to estimate the age of a system with measured 
eccentricity. Equivalently, we can express the mean induced eccentricity in a system of age 
fg 10 ® yr as 

e ~ Max I P]/^ 4 exp 



where we have defined 77 = tgnijvio. Since the cluster parameters n and v are usually 
known quite accurately, any uncertainty in 77 comes primarily from the age of the system. 
The steep dependence of e on 77 suggests that eccentricity may be an excellent dynamical 
age estimator for LMBPs in globular clusters. This is important because the timing age of a 
cluster pulsar is usually unreliable, the measured pulse period derivative being contaminated 
by the acceleration of the pulsar in the mean gravitational potential of the cluster (Phinney 
1992). 


In Eigure 2 we show the results of applying equation ( 8 ) to the nine LMBPs currently 
known in globular clusters. Eor all systems we calculate the range of expected values for 
the eccentricity corresponding to a pulsar age between 10® yr and 10^° yr. All reliably 
measured timing ages of LMBPs (both in clusters and in the Galaxy) fall inside this range. 
Unless otherwise indicated below, all data are from the compilation by Phinney (1992, 
Tables 1 and 2). We also show for reference the contours of constant 77 in the (P, e) plane, 
calculated from equation ( 8 ), as well as the theoretically predicted “birth line,” following 
the eccentricity-period relation for unperturbed systems (Phinney 1992). With very few 
exceptions Galactic LMBPs are observed to lie close to this line (Phinney & Kulkarni 1994). 


Eor six out of nine currently known LMBPs in clusters, the orbital eccentricity has not 
yet been measured. Eor these systems, our results can be taken as crude order-of-magnitude 
predictions. We note that all these systems lie in the region of the plot where the use of 
the exponential result (eq. [4]) would lead to completely incorrect results (underestimating 
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the predicted e by many orders of magnitude). Note also that all systems have a predicted 
eccentricity that is in principle measurable: values as small as e ~ 10“® are detectable by 
pulsar timing (PSR J2317+1439 has a measured eccentricity of 1.2 X 10“®; Camilo, Nice, & 
Taylor 1993). In all cases where an upper limit on e has been quoted by the observers, it is 
compatible with our predicted range. For the M53 binary, the current upper limit e < 0.01 
(Kulkarni et ah 1991) falls just inside our predicted range and implies a maximum age of 
9 X 10® yr for the pulsar. 

Three systems have a measured and rather large eccentricity: PSR B1516+02B in M5 
(e = 0.13; Wolszczan et ah 1989), PSR B1620-26 in M4 (e=0.025; Thorsett, Arzoumanian, & 
Taylor 1993), and PSR B1802-07 in NGC 6539 (e=0.21; D’Amico et ah 1993). The last two 
systems are thought to have more complicated dynamical histories than we have considered 
here. The M4 pulsar is actually in a hierarchical triple system which was probably formed 
through an exchange interaction involving the LMBP and another binary (Rasio, McMillan, 
& Hut 1995). The inner binary pulsar is likely to have its eccentricity perturbed mainly 
by the second companion in the current triple (Rasio 1994). For the NGC 6539 binary, 
the very large measured eccentricity, e = 0.21, is about two orders of magnitude above 
the maximum possible value induced by distant interactions with passing stars (obtained 
from eq. [8] when t ~ 10^° yr). In addition, this system also violates the predicted relation 
between companion mass and orbital period for a standard LMBP (Phinney 1992). Thus 
the present companion was probably acquired through a more violent interaction that left 
it in a highly eccentric orbit. The M5B binary also has a very large eccentricity, e = 0.13, 
but this value is consistent with having been induced by distant interactions with other 
stars in the cluster. Using equation (7) we predict an age t = 7 X 10® yr and a pulse period 
derivative Pp < 2 X 10“^° for the pulsar [assuming t < Pp/{2Pp)]. 

We thank Sterl Phinney and Steve Thorsett for many useful discussions, and Steve 
McMillan for first pointing out to us the fundamentally different behavior of the induced 
eccentricity in the unequal-mass case. F. A. R. is supported by a Hubble Fellowship, funded 
by NASA through Grant HF-1037.01-92A from the Space Telescope Science Institute, which 
is operated by AURA, Inc., under contract NAS5-26555. D. G. H. thanks the Institute for 
Advanced Study, Princeton, for its hospitality while work on this project was being carried 
out. 
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Fig. 1.— Final orbital eccentricity ef following an interaction with a star on a coplanar 
prograde parabolic orbit with pericenter distance Tp (given in units of the binary semi-major 
axis a). The masses of the two binary components are mi = 1.4 M 0 and m 2 = 0.2 M 0 ; the 
mass of the passing star is m 3 = IM 0 . Solid lines show the analytic results (eqs. [l]-[5]); 
dots show the results of numerical integrations. Round dots are for a system with initial 
eccentricity = 10“®, squares for = 10“®, and triangles for = 10“^. The “error bars” 
show the full extent of the dependence on the orbital phase. The dashed line shows the 
exponential decay predicted by eq. (4) alone. 


Fig. 2.— Orbital period (in days) and eccentricity of the nine currently known LMBPs in 
globular clusters. The round dots show the positions of the three systems with measured 
values of both e and P. The vertical dotted segments show the predicted eccentricity for 
a pulsar age between 10® and 10^° years (eq. [ 8 ]). For the M53 binary we have taken into 
account the current upper limit (e < 0.01). Solid lines show the variation of e with P for 
fixed T] = tgu^lvio, and are logarithmically spaced with two lines per decade of 77 . The dashed 
curve shows the theoretical “birth line,” along which unperturbed systems are expected to 
lie. 








